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We perform combined resistivity and compressibility studies of two-dimensional hole and electron 
systems which show the apparent metal-insulator transition - a crossover in the sign of dR/dT with 
changing density. No thermodynamic anomalies have been detected in the crossover region. Instead, 
despite a ten-fold difference in r s , the compressibility of both electrons and holes is well described 
by the theory of nonlinear screening of the random potential. We show that the resistivity exhibits a 
scaling behavior near the percolation threshold found from analysis of the compressibility. Notably, 
the percolation transition occurs at a much lower density than the crossover. 
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The apparent metal-insulator transition (MIT) in 
high-mobility two-dimensional systems remains a topic of 
fundamental interest and continuing debate 0. The 
anomaly of these systems is exemplified by the existence 
of a narrow range of carrier densities around n = n c 
where the slope of the temperature dependence of the 
resistance, dR/dT, changes its sign. To unravel a com- 
plex interplay between interactions and disorder in this 
phenomenon, it is essential to combine transport mea- 
surements with other experimental probes, in particular 
measurements of the thermodynamic density of states 
(also referred to as the charge compressibility 00) 
X = dn/d^,, where fi is the chemical potential. There 
have been only few measurements of x near the appar- 
ent MIT 0, , among which work on a 2D hole 
gas with large values of the Coulomb interaction param- 
eter r s = l/yVnaJ r* 5 — 16 has attracted much atten- 
tion. (Here as = 18 A is the effective Bohr radius for 
the hole mass of 0.38 mo.) In their experiments done 
at T = 0.3 - 1.3 K the authors of Ref. found that 
the inverse compressibility x~ 1 (?t.) has a minimum which 
is positioned exactly at n c . This was interpreted as a 
thermodynamic signature of an interaction-driven phase 
transition discussed in theoretical works . 

An alternative explanation of the minimum of X _1 ( 71 ) 
can be based on the nonlinear screening theory 

(nst) mm HQ 

that emphasizes the role of disor- 
der. The basic premise of the NST is that a low-density 
metal is unable to screen fluctuations of potential, so that 
depiction regions with vanishingly small local density ap- 
pear and grow as n decreases. The NST predicts that 
X &) has a minimum at n = n m (determined by disor- 
der), and a rapid upturn to positive values at n < n m . 

This theory also predicts a percolation threshold at 
n = n p where n p sa n m /3 in typical GaAs sys- 



tems |lj| . There have been suggestions, based on the con- 
ductance scaling, that the percolation transition is cl osely 
related to the change in the sign of dR/dT [II El Elf. 
(The existence of the percolative MIT in 2D GaAs struc- 
tures was proposed earlier in 

In this work we use combined compressibility and con- 
ductance measurements to shed light on the origin of the 
apparent MIT in 2D hole gases with large interactions be- 
tween the carriers - a problem widely debated over the 
last few years 0, 0, 0] . We compare the behaviour of 
the compressibility of holes with that of electrons (with 
much weaker interactions) and conclude that the appar- 
ent MIT in 2D hole gases is not related to a quantum 
phase transition. Similar to the case of 2D electrons, it is 
most likely caused by an interplay of different scattering 
mechanisms [l^, |2(], EJ. The combined measurements 
also enable us to demonstrate that there is no direct link 
between the MIT and the percolation transition either. 

Two types of system have been examined: a 2DEG 
with r s = 1 - 4 and a 2DHG with r s = 10 - 35. We find 
no relation between the position n m of the minimum in 
X^ 1 and the position n c of the MIT - the two densi- 
ties can differ by a factor of two and, furthermore, their 
ratio is sample and cooldown dependent. The x _1 ( n ) 
dependence fits the NST E3 predictions very well for 
both carrier types. We show that in some range of T 
the dependence of the conductance of both systems on 
carrier density fits the usual percolation scaling ansatz. 
The found percolation threshold agrees with the predic- 
tion from the analysis of the compressibility in terms of 
the NST. However, its density is significantly lower than 
the crossover density n c , which makes a direct, universal 
connection between the apparent MIT (the sign change 
of dR/dT) and percolation doubtful. 

Our 2DEG structures E01 and E02 contain two GaAs 
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FIG. 1: Capacitance of a double-layer 2DEG structure 
E01 measured simultaneously by two methods: Cb - field- 
penetration, C - capacitance; T = 4.2 K, / = 4 Hz. Inset: 
diagrams of capacitance measurements in (a) double- and (b) 
single-layer structures. 



quantum wells of width w = 200 A separated by a 200 A 
and 300 A-thick AlGaAs barrier, respectively, Fig. 2] The 
top-layer mobility at the highest density is 5x 10 5 cm 2 /Vs 
for E01 and 8 x 10 5 cm 2 /Vs for E02. The 2DHG samples 
H03 and H05 with mobility 4.7x 10 5 cm 2 /Vs and 5.6x 10 5 
cm 2 /Vs, respectively, arc formed in standard single-layer 
AlGaAs/GaAs heterostructures with Au gates. 

The resistance as a function of gate voltage, V g , is mea- 
sured at T = 0.03 — 10 K, in a similar way to [lij . 
and the relation between V g and n is established by 
the Hall effect. The position of n c is determined at 
the lowest temperature of the compressibility measure- 
ments, T = 0.26 K. The compressibility is found to be 
temperature independent in the range 0.26 — 1.5 K for 
2DHGs and 0.26 - 5 K for 2DEGs. The compressibil- 
ity and conductance measurements have been repeated 
several times, one after another, to confirm that there is 
no drift in the structure. The compressibility is deter- 
mined using two techniques - the "capacitance" and the 
"field-penetration" , Fig. 

Capacitance method — Here the ac voltage V g = 2 mV 
of frequency / = 1 — 100 Hz is applied to the gate 
and the 90° phase-shifted current I in the probed layer 
is measured (Fig. ^ inset), yielding the capacitance 
C = I/2irfV g . Instead of x we discuss the quan- 
tity which is a measure of inverse compressibility: the 
Thomas-Fermi screening radius d = (eeo/e 2 )x _1 0. It 
is easy to show that d is related to C by the formula 
d = e£o-4(C _1 - Cq 1 )/(1 + 7 _1 ), where A is the gated 
area and Co = eeoA/si is the geometric capacitance be- 
tween the gate and the probed 2D layer, Fig ^ bi a 
double-layer structure, where the "probed" 2D layer is 
the top quantum well, the factor 7 = (db + S2 + w)/si 
accounts for the electrostatic interaction between the two 
layers 0. Here S2 is the separation between the centers 



of the quantum wells and d b w ag/4 <C S2 + w is the 
screening radius of the bottom layer. In a single-layer 
structure S2 — > 00 and the correction 7 -1 is set to zero. 

Field-penetration method — In this technique the com- 
pressibility of the top layer is found from the capacitance 
Cb between the gate and the bottom layer. The measured 
capacitive current h is caused by the electric field pen- 
etrating through the top layer in proportion to its X ■ 
Once Cb(V g ) is found, the screening radius d is calculated 
using Eqs. (7)-(9) of Ref. 0. 

In real 2D systems there is a change in the transverse 
confinement (and the corresponding subband energy) 
with varying n. This gives a small contribution Ad su b to 
d whose sign and magnitude depend on the structure and 
the method of capacitance measurement. To compare the 
results obtained in different experimental situations, we 
subtract this contribution and discuss the quantity d* = 
d — Ad su b- The value of Ad su b is of the order of the 2D 
layer thickness and is calculated by a perturbation theory. 
For a single-layer heterostructure Ad su b = 0.46(a_e/n-) 1 ^ 3 
\22\ . For double-layer (quantum well) structures we com- 
pute A(i su b using the infinite square-well approximation, 
similarly to Ref. [2j|. This gives for the capacitance 
method Ad su b = +0.3967iu(l - 0.0544A) in the thin-well 
limit, A = nw 3 /a,B « 1. In the field-penetration tech- 
nique AcUb = -0.1033ti;(l - 0.1148A), cf. Ref.0. 

In Fig. |3 we present d* for sample E01 found from 
both the capacitance and the field-penetration methods. 
The field-penetration technique is more accurate than the 
capacitance technique at large n, because in this method 
the subtraction of the geometric term Co is obviated. 
However, at low densities around the minimum in \~ 1 
the contribution Co is less important, and one can see 
that the two methods give nearly identical results. 

Now we turn to the comparison of the results with the 
nonlinear screening theory. This theory predicts: 



d* = (o B /4) + Adex + Ad C or + Ad dis . 



(1) 



The first term in Eq. (|TJ is due to the single-particle 
density of states (kinetic energy) of the 2D carriers. 
The correction Ad ox = — (87r 3 ri) -1 / 2 comes from the 
exchange interaction. Another negative contribution 
is due to correlations between the carriers: Ad cor = 
(seo / e 2 )d 2 (nE c ) / dn 2 , where the correlation energy per 
particle E c is computed according to Ref. Q. As the 
carrier density decreases, Ad ex and Ad cor cause a change 
in the sign of d from positive to negative - the effect seen 
experimentally 0,0,0- Disorder, however, brings a pos- 
itive contribution, Ad,A\s , responsible for the upturn of 
X -1 at low densities pL3| : 
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where r\ = ns/y/rn, s is the spacer and n» is the effective 
2D concentration of dopants (an adjustable parameter, 
see below). 
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FIG. 2: d* for 2DEG structure E01 with a = 200 A, 
si = 700 A, S2 = 400 A. Solid lines are the results of 
the field-penetration method, dashed line - the capacitance 
method. Dotted line - the NST theory with n< = 1.2 x 
10 11 cm -2 . (The donor concentration found from T q (n) is 
n\ = 2.3 — 4.5 x 10 11 cm" 2 , while the maximum concen- 
tration of (uncorrelated) donors is nf = 9.0 x 10 11 cm" 2 .) 
Inset: Comparison with 2DEG structure E02 (s = 400A, 
si = 900 A, s 2 = 500 A), where m = 2.2 x 10 11 cm" 2 and 
n\ = 3 - 4.5 x 10 11 cm" 2 , while nf = 11 x 10 11 cm" 2 ). 



Equation was derived assuming that disorder is 
produced by a 5-doped layer of uncorrelated dopants with 
a two-dimensional concentration m |l3j . The samples 
in this work have three-dimensional doping. It can be 
shown that the impurities closest to the 2D layer have 
the greatest effect on x> an d therefore Eq. @ is still 
valid provided one uses an effective n,. Unfortunately, 
m cannot be determined in a simple way because not 
all impurities can be ionized and also because of existing 
correlations in their positions j2f|. Both factors reduce 
the effective n, [l ll| compared to a naive estimate based 
on the total number of impurities known from the growth 
conditions (denoted here by nf). A better estimate of Hi 
is deduced from the quantum lifetime r q . We have found 
T q (n) from the analysis of the Shubnikov-dc Haas effect 
and determined the effective concentration n\ according 
to Ref. [2(| - it turns out to be 3-5 times smaller than 
nf. 

The dotted line in Fig. |21 is the best fit by Eq. for 
structure E01, with rij found to be close to nj . Good 
agreement with Eq. where is close to n\ has also 
been obtained for E02 - the inset compares the results for 
the two 2DEG structures. The latter has a larger spacer 
and in agreement with the NST the X _1 ( n ) minimum, 

1 /2 

n m ~ 0.38 nf /s, is shifted to lower n. The position of 
the minimum in x _1 ( n )i indicated by an arrow, is found 
from the fitted curve with an accuracy of < 2%. (Inac- 
curacy of the calculation of Ad su b has an even smaller 
effect on the position of n m .) The range of densities 
where dR/dT changes its sign is indicated in Fig. [3 by 
n c . It is seen that the crossover in R(T) occurs at higher 
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FIG. 3: d* for 2DHG structures H03 and H05 with s = 
500 A and Si = 2670 A, obtained by the capacitance method. 
Dotted lines - NST theory. The dopant concentrations for 
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and nf = 5.6 x 10 11 cm 
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% n\ = 1.1 - 2.2 x 10 11 cm"' 
The minimum for H05 (rij = 



n. Inset: d* for H03 (cooldown 2), 
nf = 5.6 x 10 11 cm"' 



= 5.6 x 10 cm ) is shifted to larger 



1.5 x 10 11 cm"' and 



densities than the minimum in \ ( n )- 

Figure|21shows the results for 2DHG samples H03 (two 
cooldowns) and H05. The dotted lines are best fits to 
Eq. JTJ with parameters n, consistent with our analysis 
of T q (n) . In the second cooldown the upturn of d* occurs 
at higher densities and the obtained n, is also larger, i.e., 
the sample in this cooldown is more disordered. This 
supports the notion of correlations among the impurities 
since these are known to depend on thermal cycling |25j . 
As in the 2DEG structures, the apparent MIT region n c 
does not overlap with n m , although here n c < n m . This 
is not surprising because of the difference in the contri- 
butions to R(T) of electrons and holes that determine 
the position of n c . The same conclusion that n c < n m 
was drawn for another studied 2DHG sample (H06, not 
shown) . 

To understand the relation between the apparent MIT 
and percolation in 2DHGs, we have attempted to extract 
the percolation threshold n° from the fit of the conduc- 
tance to o~(n) — a(n — ji^)* following the procedure in 
13 El IH • Figure 0] shows the results with the expo- 
nent t = 2.1 ± 0.1 at T = 0.26 K for both electrons and 
holes, which is close to that in works [lol Il8j|. For the 
2DHGs the exponent decreases with decreasing temper- 
ature to t = 1.6 ± 0.1, which is close to t = 1.31 ex- 
pected for classical percolation j2||. Notably, the value 
of the percolation threshold is significantly lower than n c : 
n c jn a v = 2.9 for H03 and n c jn a v = 2.1 for H05. There- 
fore, we surmise that the apparent MIT is not due to 
percolation but, similar to 2DEGs i s due to an in- 
terplay of the "metallic" T-dependence caused by phonon 
scattering and Fermi-liquid corrections to impurity scat- 
tering, and the "insulating" dependence caused by local- 
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FIG. 4: (a) Fits of the conductance to <r(n) = a(n — rip)', (b) 
An example of the temperature dependence of the resistivity 
of H03 structure. The range of n c in this cooldown is indicated 
by the dotted line. 
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We want to stress, however, that using the percolation 
scaling in a quantum system requires justification. One 
can rationalize this procedure by arguing that the perco- 
lation concept can apply in some intermediate range of 
T, high enough to preempt localization at n > rip yet 
low enough to inhibit thermal activation at n < rip. An- 
other point of concern is that the scaling is observed over 
a broad range of densities (up to n/n° ~ 2, similar to 
|14lll5|) 1 w hile it is expected to work only near the critical 
point |27j . However, our combined transport and com- 
pressibility measurements allow us to examine directly 
the applicability of the scaling procedure and demon- 
strate its validity. Using the NST prediction [lit Il3j . 
we obtain independently n p w 0.12y / n7/s w n m j?> and 
compare it with rip. We have established that the per- 
colation thresholds derived from both methods are very 
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H03, and n p = 1.2 x 10 1 " cm" 2 and n a p = 1 x 10 lo cm 
for H05. (This comparison is done at T = 0.26 K and 
even better agreement is obtained at T = 30 mK as rip 
increases with decreasing T by about 20%.) 

In summary, our combined conductance and compress- 
ibility measurements suggest that the apparent MIT in 
2DHGs with r s up to ~ 30 is neither an interaction- 
driven phase transition nor a percolation transition. The 
behavior of the compressibility at low hole densities is 
well described by the nonlinear screening theory. This 
indicates that the upturn in x _1 (n) is due to depletion 
regions in the channel, with total area less than 3.5%, 



caused by disorder. 
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